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Abstract: Determining the spin of new particles is an important tool for discrim- 
inating models beyond the Standard Model. We show that in case of cascades of 
subsequent two body decays the existing strategy to extract the spin from lepton 
and quark spectra can be used without changes even if one aUows for dim-5 and 
dim-6 operators which might be induced by physics just beyond the reach of LHC. 
We show analytically that these operators do not change the overall structure of 
these spectra. 
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1. Introduction 

The Large Hadron Collider (LHC) has started the direct exploration of physics at the 
TeV scale and is searching for new particles which are predicted in various extensions 
of the Standard Model (SM). Many of these models predict partners of the known SM 
particles, which usually have the same quantum numbers and properties but for the 
mass and the spin. For example, in supersymmetric models the fermions have scalar 
partners whereas in models with extra dimensions fermionic partners are predicted. 
Therefore, immediately after the discovery of new particles the question will arise 
how to discriminate between the various models and spin determination will play a 
crucial role here. 

The question on how to determine the spin of new particles has been addressed 
by several authors: in [1, 2, 3, 4, 5] s-channel resonances have been investigated 
and in [6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19] spectra of cascade decays of 
subsequent two-body decays have been used to obtain information on the spin. An 
additional possibility to get information on the spin is cross section measurements 
provided one knows the representation of the particle produced [20], e.g. whether it 
is a colour triplet or a colour octet. Also three-body decays have been investigated 
in this direction for specific scenarios [21]. In ref. [22] a strategy has been worked 
out for scenarios where three-body decays are dominating. 

For decay chains of at least three subsequent two-body decays involving with 
at least three visible SM particles in the final state one can construct sufficient 
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many kinematical variables to determine the masses of four unknown particles even 
if the lightest of them escapes detection, see e.g. [18] and references therein. This 
is in particular important if one requires that a SM extension explains the observed 
dark matter (DM) due to heavy weakly interacting particles, e.g. the neutralino in 
supersymmetric models or the lightest KK-excitation of the photon in models with 
extra dimensions. However, these kinematical variables do not give any information 
on the spin of the particles involved whereas the angular distribution of the decay 
products is highly effected by the spin of the intermediate particle and each spin has 
its specific decay signature. This property has widely been studied [23, 8, 9, 10]. 
Potential complications in this context are at the one hand that one does not know 
a priori which particles belong to which decay chain. On the other hand it is usually 
also not known in which particular decay chain a certain particle has been produced 
which can potentially wash out effects as one has to sum over different possibihties. 
An example of the later case is in case of supersymmetric models the chain 



where g is a scalar quark, / is a scalar lepton and Xi 2 neutralinos. Here one can 
build kinematical variables using the jet originating from outgoing q and the two 
leptons [24]. 

The fact, that LHC has failed so far to detect physics beyond the SM might 
imply that part of the corresponding spectrum is somewhat beyond the reach of 
LHC, e.g. gluinos or KK-excitations in the multi-TeV range. They might however 
induce dim-5 or dim-6 operators which are only mildly suppressed which potentially 
endangers the above strategies to determine the spin of new particles in cascade 
decays. However, as we will demonstrate below, this is not the case and one can use 
the strategies developed so far with changes even if these additional operators are 
presented. In the next section we will fix the notation and discuss the potentially 
dangerous operators. In Section 4 we derive the main results and then conclude in 
Section 5. 

2. Setup and Kinematics 

These methods rely on the fact that one can decompose the matrix element of a 
scattering process a + b^i^c + d such that the dependence of the scattering 
angle 9bc is solely expressed in Wigner-d functions and a remaining reduced matrix 
element which only depends on the total angular momentum J, the helicities of 
the involved particles, and the particle masses. Using crossing-symmetry, one can 
rewrite the decomposition of this matrix element [25, 26] for the s-channel process 
a + b^i^c + d into the matrix element for the decay a ^ b + c + d which reads 
in the rest frame of the intermediate particle i [27] 



,0 
1 




XiXf 



(2.1) 
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with the differences of the initial and final particles' helicities Aj = (A^ — \b),^f = 
(Ac — Xd) and where is the scattering angle of two of the products. The Wigner-d 
rotation matrices 'iAi,A/ depend polynomially on cos6'bc/2 with degree 2s [28]. After 
squaring, this leads to a polynomial in cos^^c with degree 2s which gives us the 
spin-dependence of the angular distribution for this decay. Since we want to study 
frame- independent variables, we rewrite eq. (2.1) in terms of invariant mass Sff of two 
visible decay products and the cos 9hc dependence of the matrix element translates 
into 



2s 



d r{X ^ ffY) = d SffY^bi s}f (2.2) 

i 

for the differential decay rate with maximal degree 2s where we have used the fact 

Sff<xEiE2-\\p\\q\coseBC (2.3) 

with p, q, El, E2 being momenta and energies of the decay products. 

We use use following notation for the masses and momenta in the investigated 
decays 

X(p,m^) f(qi,mf) + I(pi,mi) f(qi,mf) + f(q2,mf)+Y(q3,my). (2.4) 

where X, I, Y can either be scalars, vectors or fermions. We use the usual definition 
of the Mandelstam variables 

«// s = {Qi + (12^ = {p- qzf] = (gi + q^f = (p - ^2)^; 

i = (g2 + ?3)' = (p-gi)'=:P7- (2.5) 

where s// is the invariant mass of the two visible SM fermions and is the variable 
we are interested in. If one wants to derive the differential decay rate dV/d Sff one 
usually replaces one of the variables, e.g. u via Mandelstam's relation 

u = 2 mj + ml + rriy — t — Sff (2.6) 

and integrates out the remaining invisible invariant mass, in this case t where the 
upper and lower bounds are 

t± = {m} + ml + ml-Sff) (2.7) 

\ Sff {sff - 4:mj) (mi + {ml - Sfff - 2ml (^^ + *//)) 
± ^ ^ ■ ^2.8) 

The remaining invariant mass Sff has then the kinematical bounds 

Sff_ = (2m/)2; Sff^ = (m^ - myf. (2.9) 
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Our interest are subsequent two body decays keeping track of the polarisation 
information which is transferred via the intermediate particle /. For this we use the 
narrow width approximation (NWA) as worked out in ref. [29]. The endpoints of 
the invariant mass distribution s// after using NWA get changed due to the on-shell 
condition for the particle / and reads 



/o \2 {ml - mj){mj - ml) 
= (2m/) Sff^ = -2 . (2.10) 

fit J 

3. Dimension 5 and dimension 6 operators 

In this work we are interested in the question if dim-5 and dim-6 operators can 
invalidate the spin analysis for new particles, e.g. by increasing the highest power of 
Sff. This is only possible if there are momentum dependent interactions and, thus, 
we will restrict ourselves to this subset. Moreover, we consider final states containing 
SM fermions giving further restrictions. 

The basic operator structures are given in [31] which, however, involves only the 
SM fields as external particles. In our examples new particles are allowed including 
additional gauge interactions with covariant derivatives — — lA^j^ and field 
strength F^^, — {d^A^, — di^A^ + ...). We assume that these additional gauge groups 
are broken at a scale as a new scalar gets a vacuum expectation value v^j) and 
inducing a dim-5 from a dim-6 operator via 

Having this in mind we get two classes of operators: fermion-fermion-vector (f-f-V) 
and fermion-fermion-scalar (f-f-S) . For the (f-f-V) interactions we have 



>Cf)4 =A^V'i7'' {giPi + 9rPr) ip2 + h.c. (3.2a) 
1 



J^m =i- (^1^'^" M + O'rPr) i^^F^. + h.c.) (3.2b) 



>Cd6 =^ {{DM (blPl + brPr) ^2Dy + /l.C.) 

^^m = ^^ (Wh) {biPi + KPr) M-^^n + h.c.') (3.2c) 

JCd6 =J^{Tl (ciPl + CrPr) D^lP2D^(f) + h.C.) 

^^D5 = j^ (V^ idPl + CrPr) 0^'^^') + h.c) (3.2d) 

and for the (f-f-S) interactions 

^{<t>ih {niPi + rirPr) ^2 + h.C.) (3.3a) 
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>Cd5 = ^ (0D;A)7'^ {yiPi + y^n) ^^20 + /i.e.) (3.3c) 

>C|)5 = (^7^ (-^^^^ + -Sr-Pr) (^M^2)</' + h.c) (3.3d) 

Obviously not all of them are independent, e.g. by partial integration one can transfer 
one of the derivatives to the other two fields. In principle one could also consider dim- 
6 operators for the (f-f-S) case but as we will discuss below, no additional features 
will occur in such a case. 



4. Results 

We will first the discuss the impact of operators containing two fermions and scalar 
as here the analytical formulas are rather compact. Then we will turn to decays 
involving also vector bosons. 

4.1 A simple example 

We start with a simple example, namely the case where X and Y are fermions and 
/ a scalar and we denote this case for later use by {FF)s. In this case that there 
are no spin correlations between the SM fermions, the matrix element, including the 
anomalous couplings due to the dim-5 operators, is of the form 

Md4,d5 oc (^(/)(1 + XiI^M^)) + A2VX/)) (4.1) 

where we have used a combinations of and £|)5, e.g. the last two operators in 
eq. (3.3). The ki are weighted sum of the fermion momenta, e.g. ki = aip + (5iqi and 
k2 = Oi2q2 + /9253, and Aj are measures of the relative importance of the dim-5 with 
respect to the dim-4 operators. Here we have put for simplicity scalar couplings only. 
As all fermions are on-shell, this immediately implies that one can use the equation 
of motion and replace in the }f,i the momenta by the corresponding fermion masses. 
This immediately implies that no additional power of s// occurs and we get 

4^ = ^ . 5° = A. (4.2) 

In the case that is involved one has to use the momentum conservation of 
the vertices to express the momentum of the scalar by combinations of the fermion 
momenta. One can also show easily in this case that using of the equations of motion 
gives only additional mass terms. 

The more comphcated case is that X and Y are scalars and / a fermion as in 
the case / relates the polarisation of the two SM fermions and we denote this kind of 
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j 



Figure 1: Feynman diagrams for the decays involving intermediate fermions and a 
scalar/vector boson. The boxes denote combinations of dim-4 and dim-5 couphngs. 



decay for later use as (SS). However, in this case also an explicit calculation shows 
that no additional power of s// occurs and we get 

dr 



ds 



ff 



A + B-Sff 



(4.3) 



The result for the coefficients A and B is given in the appendix. One can show along 
the same lines, that operators of the form 



(4.4) 



also do not give any contributions to these two decay chains. As we only wanted to 
indicate the principle structure of the couplings we did not write any chiral couplings 
as they are not important in this context. 



4.2 Decays involving vector bosons 

Here we have four cases, three where the intermediate particle is a fermion 

S ^ ffV 
V ^ ffS 
V^ffV 

yielding a differential decay rate of the form 

d r 



{SV) 
(VS) 
(VV) 



ds 



ff 



(4.5a) 
(4.5b) 
(4.5c) 



(4.6) 



where A, B contain the dependence on masses and couplings involved. The Feynman 
graphs for the corresponding decays are the second to fourth in Fig. 1. The final 
decay process shown in Fig. 1 contains an intermediate vector boson 



V 



ffF 



{FF) 



V 



yielding a partial width of the form 

dr 



ds 



= A + B ■ s}f + C ■ s"^ 



ff 



ff 



(4.7) 



(4.8) 
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A, B, C are again functions of the involved couplings and the masses. In all cases we 
have already anticipated that no additional powers of sjj are induced. The only effect 
of the anomalous couplings is to change somewhat the relative size of the coefficients 
A, B and C. 

We sketch the main results for the case of {VV)f as this is the one where one 
might expect the highest power of momenta and thus the highest power of s//. The 
amplitude for this decay leads to the following matrix element 



(4.9) 



where we have indicated the parts stemming from dim-4 (dim-5 ) operators by T^y 
(rf y), where i — 1,2 indicate the first or second vertex in the decay. For the squared 
matrix element we get the following structure 

pols pols 



ml ml 



~ -4 [ml - 2 p-qiY {ml qi ■ q2 - 2 p ■ qi p ■ q2) (4.10) 

using only dim-4 contributions and where we have omitted all chiral couplings as they 
are not important at this stage. After taking the NWA, the Mandelstam variables 
become 

* = (p - gi)' = (^2 + gs)' = (4.11a) 
{qi + qzf = (p - q2f = 2mj + ml + m^ - Sff - mj (4.11b) 

Therefore, in NWA we can only get a s// term from the scalar products 

p ■ q2 ^ Sff + const.; qi • q^ ^ Sff + const. (4-12) 

while the remaining scalar products just give a mass squared and hence constant 
contributions. Therefore we get as claimed a linear function in s// after taking the 
NWA. Similarly one finds the s// dependence in the cases {SV) and {VS) including 
only the dim-4 operators. It turns out that the impact of the dim-5 operators on 
these decays is the same as for the {FF)y case and, thus, we will present first the 
main facts for this decay using dim-4 operators only and discuss the dim-5 operators 
for all four decays taking the (VV) case as an example. 
In the case of {FF)v, the most involved one, we get 

= i^?v + ^fv),u{p)) ■ {u{qs) {r^Ur^^)Aq2)) 
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indicating again the parts coming from dim-4 (dim-5 ) operators by Tfy (rfy). For 
the squared amphtude we obtain for the part stemming from the dim-4 operators 

pols 
pols 

mj mj 

~ IGi-g'^'^'p ■ q, + p>^q^' + p%'){-g'^^' q, • ^3 + q^^q^^' + q^q^^) 

xir-^^W'''-^^) (4.14) 

In comparison to the case (VV) we have the additional contributions of p^q^ . . . 
which can then be contracted with the momenta combination of the trace of the 
second fermion hne. These contractions give new momenta combinations and one 
can see that terms of the form (^i •52) {P'qs) arise which give according to eq. (4.12) a 
0{s'jj-) contribution which is the characteristic for intermediate vector bosons. Note, 
that {p ■ qs) = Sff + const, according to eq. (2.5). 

In both cases the couplings induced by the dim-5 operator have in the most 
general case the form 

rf^'' = Ai ■ a^^kv a + Bi k't, + Q ki^, (^ = 1, 2) (4.15) 

where ky is the momentum of the vector boson, and ku {k2i) are linear combinations 
of the particle momenta at the first (second) vertex. Plugging this in eq. (4.9) and 
considering the part of the matrix element squared with the anomalous coupling 
squared, as this gives potentially the highest power in Sff, we find 



|2 



^ \MvV, D5\ 
pols 

J2 {u{qi){Aia''''pa + B,p>i + C^q^MMa'^'q^ a + ^sp^ + C^q^Hq^)) 



pols 



X (^Z(5i)(^ia^'>„ + B^p>^' + Ciq'^yjiA2a^\3 a + B^p'f + C2qi)v{q2))^ 
X eM(p)v(p)*e.(g3)e.'(93)* • (4.16) 



^ V- 
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We see immediately that the a^^'^ terms drop out as they are contracted with sym- 
metric products of momenta. The remaining products of momenta are 



PiPn = (P - QiTPh , QiPn , p'}Q3,h 



{q2 + q3Tq3,n , g2?3„ 



(4.17) 



Using eq. (4.11) wc sec that in NWA all of them give only sums of masses squared 
but no additional powers of Sfj. The same reasoning can also be used to show that 
also in the case of (SV), (VS) and {FF)v no additional powers of s// are induced. 

The main reasons, why higher dimensional operators do not change the overall 
lepton and quark spectra of the decays, can be summarized as follows: 

1. Additional S// dependence, which is equivalent to additional cos^^ dependence, 
can only arise through the following products: {qi ■ q^) or {p ■ 52)- AH other give 
in NWA sums of masses squared, see eqs. (4.11a) and (4.11b). 

2. The antisymmetric part, e.g. the a^"^ part, of the (f-f-V) coupling gets always 
contracted by the same momentum due to the polarization sum/propagator of 
the vector boson and hence gives zero. 

3. The momentum dependent parts in the (f-f-V) coupling relate only momenta 
within a given vertex. In NWA the momentum conservation at a given vertex 
implies that all scalar products of momenta can be expressed either as masses 
squared or as t = mj. 

4. Momenta contracted with gamma-matrices yield only masses after using the 
Dirac equation. 

Therefore, also dim-5 operators where fermions are coupled to vector bosons do not 
change the highest power in Sff. 

We have checked that the same reasoning also applies for dim- 6 operators. Higher 
than dim-6 operators should play no role as latest at this stage higher order correc- 
tions due to emission of gluons and photons become more important. 

4.3 Numerical example 

We have checked numerically several examples to test the quality of our reasoning 
in the previous sections. As a random example, we show for all six decays the 
differential decay rate f for the dim-4 and dim-5 result in fig. 2, The masses and 
the couplings are chosen as 



= 1 TeV 
m/ = 0.4 TeV 
ni — IjUr — 0.1; 
ai — l,ar — 0.3; 



gi = l,gr^ = 0.1; 
bi = 0.4, br = 0.5; 



ruy = 0.15 TeV 
A = 3 TeV 



ci = 0.8, Or = 0.7; 



ruf — rrif 



= 



(4.18) 
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D4W/D5 






D4 w/o D5 - - - - 






D4 w/ D6, NWA 






D4 w/o D5, NWA 
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S„/TeV^ 







^^S) 




















D4W/D5 




D4 w/o D5 - - - - """VV^. 




D4 w/ D5, NWA '"Xl 




D4 w/o D5, NWA \ 
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D4 w/o D5 - - - - 






D4 w/ D5, NWA 






D4 w/o D5, NWA 
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D4 w/o D5 - - - - 






D4 w/ D5, NWA 






D4 w/o D5, NWA 





0.3 0.4 

s,/TeV^ 




D4W/D5 ■ 
D4 w/o D5 ■ 
D4 w/ D5, NWA ■ 
D4 w/o D5, NWA ■ 
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S„/TeV^ 



(FF)s 



D4W/D5 ■ 
D4 w/o D5 ■ 
D4 w/ D5, NWA ■ 
D4 w/o D5, NWA ■ 



0.3 



0.4 



Figure 2: Differential decay rate over the invariant mass sjf for dim-4 couplings 
only (black) and including additional dim-5 (red/grey) couplings for intermediate 
fermions/bosons (first two row/last row). The lines correspond to: full red lines: ex- 
act spectra including dim-4 and dim-5 operators, dashed lines: exact spectra taking dim-4 
operators only, full grey lines: spectra including dim-4 and dim-5 operators using NWA, 
dotted lines: exact spectra taking dim-4 operators only using NWA. 

xi = 0.4, Xr = 0.3; yi = 0.5, y,. = 0.3; zi = 1, Zr = 0.1 . 

We show for each decay the complete dim-4 (black dashed) and dim-5 distribu- 
tion (red sohd) and for comparison the NWA results (dim-4 , NWA: black dotted; 
dim-5 , NWA: grey solid). For all decays, the NWA result is a very good approxima- 
tion for the complete decays including off-shell effects. The main differences are at the 
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kinematical endpoints. The influence of the dim-5 operators in this mass/coup hng 
scenario is the largest in the {FF)y followed by the {VS) decay. The flrst one is 
because in this case we have the highest power in s//, see eq. (2.2), and is a quite 
generic feature. The second case depends much stronger on the couplings and masses 
of the scenarios considered. For completeness we note, that for {FF)s no changes 
are visible when comparing the cases with and without dim-5 operators, because we 
show the normalized distributions. 

5. Discussion and conclusions 

The determination of spins of new particles is an important task at the LHC once 
an extension of the SM is discovered. In this way one can also discriminate between 
model classes, e.g. between supersymmetry and extra dimensions. In case that there 
are two subsequent two-body decays involving two SM fermions one can show that 
the spin of the intermediate particle reflects itself as in the highest power of s// in 
the differential decay width dV/dsff. 

There are two main theoretical uncertainties in this kind of considerations: real 
emission of photons and gluons in the decays. These can be calculated once the 
quantum numbers of the new particles are know. The second class are higher than 
dim-4 operators which are induced by particles close to the LHC reach. In this 
paper we have shown that such operators do not induce higher powers of Sjf in 
the partial width and thus the existing analyses are valid even if such operators are 
numerically important. However, their presence can change the exact form of the 
slope considerably as we have seen in a concrete example. This potentially affects 
the determination of the underlying parameters. 
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Appendix 

As an example for the NWA result of the decays including anomalous couphngs, we 
give here the decay width for S — > ffS. We used the following short forms for the 
couplings from eq. 3.3 and the masses used 



Tl = 



mi 




my 




mx 



rrix 
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9v^ 9r + Qh QA^Qr- Qh nv^nr + nf, UA^Ur- ni; 

and analogous for ai/j.- bi/r, ci/r, xi/r, yi/r, zi/r- The NWA result in terms of these is: 

111 



with 



d Sff {27r)^ 32m X rniT I 



Ada = 2m^ (l - rf) {rf - Ty) {{g^ + Qy) {n\ + ^v) + "^dAgvnAnv) ] 
Bd4 = - ISfi-A gv nA riv r|; 

Ad5 = - 2m^r| {rf - l) {rf - 4) x [a\R',^ (nl + 

+t| (^RJ {x\ + X^) - 2R'^R'^{xaZa + XyZy) + {z\ + A) 

+2aAR'a {^(^{civK - cvR'JiKxA - R',za){Kxv - R'.Zv) 

-CaR'c (^R'^ {xa + Xy) - 2R'^R'^{xaZa + XyZy) + RJ {z\ + 

+2avnAnvR'a - R'c {^A {n\ + riy) + 2cynAny) j 

+r] [a^R'^ (<2 (^^ + 4) _ 2R'^R'^{xAZA + xyZy) + RJ (4 + 4)) 

- 2avR'aR'c[2cA{KxA - R',za){Kxv - RW) 

+Cv [R^^ (4 + Xy) - 2R'^R'^{xaZa + XyZy) + RJ (4 + 4))) 

+RJ ({c\ + Cvf (^R'x (4 + 4) - '2R'xK{xAZA + XyZy) + RJ (4 + 4)) 

+ 4:CACy{R'^XA - R'^ZA){R'^Xy - R'^Zy)^ 

2 2 '2 2 2 '2 / / 

+ ayUARa + ayUyR^ — AayCAUAnyR^R^ 

- 2ayCynlR'^R'^ - 2ayCynlR!^R!^ 

+ 4^A-R? + 4^y-R? + AcACyUAnyR^^ + Cyn\R'^ + CyUyR!^ 

+ (4 + 4) {r'x {xl + 4) - ^R'^KixAZA + XyZy) + RJ {z\ + 4)) 

+ AgAgviR'xXA - R'^za){R'xXv - Kzv) ; 
Bd5 = - 164 {{(^aK - CaR'c) {avK - cvK) {^Any 

+Tj{R',^Xa — R'^ZA){R'xXy — R'zZy)) 
+gAgv{R'xXA - R'^ZA){R'xXy - R'^Zy)) . 
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